Correlation Assisted Phonon Softenings and the Mott-Peierls Transition in VO2 



Sooran Kim, Kyoo Kim, Chang- Jong Kang, and B. 1. Min 
Department of physics, PCTP, Pohang University of Science and Technology, Pohang, 790-784, Korea 

To explore the driving mechanisms of the metal-insulator transition (MIT) and the structural 
transition in VO2, we have investigated phonon dispersions of rutile VO2 {R-YO2) in the DFT and 
the DFT+[/ {U: Coulomb correlation) band calculations. We have found that the phonon softening 
instabilities occur in both cases, but the softened phonon mode only in the DFT+?7 describes 
properly both the MIT and the structural transition from R-YO2 to monoclinic VO2 (M1-VO2). 
This feature demonstrates that the Coulomb correlation effect plays an essential role of assisting 
the Peierls transition in i?-V02. We have also found from the phonon dispersion of M1-VO2 that 
Ml structure becomes unstable under high pressure. We have predicted a new phase of VO2 at high 
pressure that has a monoclinic CaCb-type structure with metallic nature. 

PACS numbers: 63.20.dk, 63.20.D-, 71.30.+h, 71.10.Fd 



Vanadium dioxide (VO2) is one of the most explored 
transition metal oxides due to its intriguing metal- 
insulator transition (MIT) and the concomitant struc- 
tural transition upon cooling. At ambient pressure 
and high temperature, VO2 has a tetragonal rutile-type 
structure {R-VO2) with metallic nature. Upon cooling 
below 340 K, R-VO2 undergoes the structural transi- 
tion to a monoclinic structure (M1-VO2) with insulating 
nature. [1, 2] The mechanism of MIT in VO2 has been 
a longstanding subject of controversy. In the structural 
transition from R-VO2 to M1-VO2, V ions construct the 
dimerization and the zigzag distortion. In i?-V02, V ions 
are centered at the distorted Oq octahedra, which are 
edge-shared along the c-axis. Due to the crystal field, V 
3d states are split into aig (g?||), (tt*), and states 
in order of energy, and so one electron of ion oc- 
cupies the lowest aig state. While the zigzag distortion 
increases the energy of bands, the dimerization of V-V 
causes the splitting of aig bands to open the gap at the 
Fermi level (£^f)-13| This kind of structural distortion is 
explained by a typical Peierls transition. 

However, the density functional theory (DFT) band 
approach fails to describe the insulating nature of Mi- 
VO2 properly. The energy gap at Ef can be ob- 
tained only when the extra Coulomb correlation U ef- 
fect of V 3d electrons is considered, which indicates that 
M1-VO2 is a Mott-Hubbard type insulator. Hence the 
Mott-Hubbard transition was proposed as the mecha- 
nism of MIT in VO^Ji-Ii Beyond the DFT band ap- 
proach, the GW(9Mll| or the hybrid functional band 
method [12] was employed to describe the insulating na- 
ture of M1-VO2 properly. There were also the DMFT 
(dynamical mean-field theory) approaches to explain the 
insulating nature of M1-VO2 bv considering the Mott- 
Hubbard U explicitly. [5|, |13H15| Therefore the consen- 
sus at the moment is that some amount of U is neces- 
sary to describe the insulating nature of Mi -V09.[16l.[l7t 
However, whether R-VO2 is a strongly correlated sys- 
tem or whether U is necessary for the MIT has not 
been clarified yet, despite a few existing studies using 



the DFT+[/pJ,|I9[ and the DMFT.[5|, UJ, |20( 

Besides i?-V02 and M1-VO2, other structural types 
of VO2 were reported, such as monoclinic M2-VO2, M3- 
VO2 and triclinic r-V02, which are stabilized under the 
uniaxial stress or with doping of Cr or Al. 21-25| All the 
monoclinic phases of Mi, M2, and M3-VO2 are insulators 
at ambient pressure, but at the pressure above 10 GPa, 
M1-VO2, Cr-doped M2-VO2, and Cr-doped M3-VO2 be- 
come metallic with some type of monoclinic structure, 
so called, Mx phase. 2^2^ Still the real structure of Mx 
phase has not been identified yet. 

As described above, there have been extensive elec- 
tronic structure studies on VO2. By contrast, there have 
been only several phonon studies on VO9 . [29l43ll. [33I. Isit 
Especially, there has been neither experimental nor the- 
oretical report on the phonon dispersion curve for VO2 
yet. Since the Peierls transition is closely related to the 
phonon softening instability, the study of phonon disper- 
sions of VO2 is expected to give a clue to understanding 
the mechanism of MIT in V02.[32|-[33 

In this letter, we have revisited the MIT and the struc- 
tural transition of VO2 by investigating the phonon dis- 
persions of R-VO2 and M1-VO2. We have found that 
i?-V02 is a strongly correlated system with U > 4.0 eV, 
and the Coulomb correlation effect plays an essential role 
in the structural transition from R-VO2 to M1-VO2. We 
have also studied the structural stability of M1-VO2 un- 
der pressure, and found that Mi phase becomes unstable 
to a phase that seems to be related to Mx phase. 

We consider the two phase of VO2, R-VO2 and Mi- 
VO2, the space groups of which correspond to PA2/mnm 
and P2i/c, respectively. [3] For the electronic structure 
and phonon dispersion calculations, the pseudo-potential 
band method and the supercell approach that are im- 
plemented in VASP[37] and PHONOPY[38] are used. 
The utilized exchange-correlation functional is the gen- 
eralized gradient approximation (GGA). The adopted 
values of U and J in the DFT+/7 are 4.2 eV and 
0.8 eV, respectively. Q All the phonon calculations were 
done after the full relaxation of the volume and atomic 
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FIG. 1: (Color online) (a) The phonon dispersion curves of 
i?-V02 in the DFT. The inset figure shows the normal mode 
of softened phonon at F. (b) The phonon dispersion curves 
of R-YO2 in the DFT+[/. Both figures show the negative 
phonon frequencies (imaginary frequencies) corresponding to 
the phonon softening instability. 



positions. [39] The initial structural parameters of R- 
VO2 and M1-VO2 before the relaxation are taken from 
McWhan et. a/.[3j] and Longo et. a/.[40] respectively. 



i?-V02 is stable only at high temperature. The most 
noteworthy is the marked difference in the phonon dis- 
persion curves between the DFT and the DFT+/7. In 
the DFT, the phonon softenings occur at q = F, M and 
X (Fig. ma)), while, in the DFT+/7, they occur at q = 
R, A and Z (Fig. Hl^b)). Indeed, the phonon softening 
at q = R was once predicted to be responsible for the 
transformation from R-VO2 to M1-VO2. 32-36| 



We have examined the normal modes of the softened 
phonons. The normal mode at q = F in the DFT is plot- 
ted in the inset of Fig. [U^a), and those at q = R in the 
DFT+/7 are plotted in Fig.E^a) and (b). Recall that the 
main structural changes from i?-V02 to M1-VO2 are the 
dimerization and zig-zag distortion of V ions, as shown in 
Fig.[2fd). It is evident that those lattice distortions can- 
not be described by the normal mode at F in the DFT. 
The softened modes at M and X in the DFT do not de- 
scribe the structural distortions either. In contrast, the 
normal modes at R in the DFT+/7 are consistent with 
the lattice distortions of VO2. As shown in Fig.[2fa) and 
(b), there are two degenerate softened phonon modes at 
R. The first mode in Fig. [2fa) represents the dimeriza- 
tions of half of V ions and the orthogonal displacements 
of the other half. The second mode in Fig.[2fb) is just the 
reverse of the first one. Note that the mode predicted by 
Gervais et a/.[36[ is close to the first mode. Also similar 
modes to above two were once obtained by using a simple 
interatomic potential model. [42] A linearly superposed 
mode using these two normal modes in Fig. [2fc) reveals 



We have first obtained the electronic structures of R- 
VO2 in the DFT and the DFT+[/. In the DFT, stable 
metallic phases are obtained both in the nonmagnetic 
and the magnetic band structure calculation, in agree- 
ment with the experiment. However, in the DFT+t/, 
more stable insulating phase is obtained in the mag- 
netic band structure calculation, [41] which is seemingly 
in disagreement with the experiment. In fact, the insu- 
lating phase of R-VO2 in the DFT+/7 had been reported 
before. [18., ilQj This discrepancy can be resolved by con- 
sidering the competition between the magnetic instabil- 
ity and the structural instability in R-V02- Of course, in 
nature, the structural instability wins over the magnetic 
instability, and so R-VO2 undergoes the structural tran- 
sition to M1-VO2 rather than the magnetic transition. 

Figure [1] shows the phonon dispersion curves of i?-V02 
in the DFT and DFT+/7. Four Raman modes (Bi^, E^, 
Ai^, and B2^) are obtained in both cases, in agreement 
with experiments. 29|-3lj The phonon softening instabil- 
ities are obtained both in the DFT (Fig. [Il^a)) and in 
the DFT+/7 (Fig.H^b)), which imply that R-VO2 is not 
a stable structure. These results are expected because 
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FIG. 2: (Color online) The normal modes of the softened 
phonons of i?-V02 in the DFT+[/ at q = R. There are two 
degenerate modes, as shown in (a) and (b). (c) A linearly su- 
perposed mode at q = R using two degenerate modes of (a) 
and (b). The blue and red balls represent V and O, respec- 
tively, (d) Schematic picture of V-V pairing and distortion in 
M1-VO2. Blue, green, and black dots correspond to V ions in 
M1-VO2, R-YO2, and M2-VO2), respectively. 
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FIG. 3: (Color online) (a),(b) The Fermi surface of R-YO2 
in the DFT and in the DFT+[/, respectively. (c),(d) x(q) of 
R-VO2 in the DFT and in the DFT+[/, respectively. (e),(f) 
Local charge density near Ef around V ion (integrated over ~ 
-1.2 eV to Ef) in the DFT and in the DFT+[/, respectively, 
in unit of (e/A^). 

simultaneous dimerizations and zig-zag distortions of V 
ions, in good agreement with the observed lattice dis- 
tortions in M1-VO2 (Fig. EJd)). The other superposed 
mode also reveals similar displacements of V ions. This 
agreement demonstrates that the softened mode at R in 
the DFT+/7 describes the structural transition of VO2 
properly. The softened modes at A and Z in the DFT+/7 
are also related to the dimerizations of V ions, but not 
directly to the zig-zag distortions of V ions. 



The role of the Coulomb correlation effect in the struc- 
tural transition in VO2 has not been invoked as an es- 
sential factor before. The present phonon study clearly 
demonstrates that the Coulomb correlation in R-YO2 fa- 
cilitates the Peierls-type structural transition. Namely, 
in i?-V02, the Coulomb correlation effect and the Peierls 
distortion are mutually cooperating in driving the MIT 
and the structural transition. To examine the U effect 
in more detail, we have checked the phonon dispersion 
curves with varying U. With increasing [/, we have found 
that the softenings at q = T, M and X disappear, whereas 
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FIG. 4: (Color online) The phonon dispersion curve of Mi- 
VO2 at the pressure of 35.1 GPa in the DFT. The inset shows 
the normal mode of the softened phonon at q = F. The blue 
and red balls represent V and O, respectively. 

the softenings at q = A, R and Z emerge for U > 4.0 eV. 
This feature indicates that large enough U is required to 
produce the right Peierls distortions in VO2. 



Figure [3] shows the Fermi surface, the electric suscepti- 
bility x{^)i arid the local charge density of R-VO2 in the 
DFT and the DFT+/7 calculations. By comparing the 
Fermi surfaces in Fig. [3^a) and (b), one can notice that 
the Fermi surface becomes flatter in the DFT-ht/, sug- 
gesting a possible nesting effect along the c-axis.ji^ This 
change is more apparent in x(q)- As compared to x(q) 
in the DFT, x(q) in the DFT+/7 exhibits high values at 
q=(x,0,0.5) {0 < X < 0.5), which explains the origin of 
phonon softenings at q=(x,0,0.5), that is, at Z and R in 
Fig. [11(b). Local charge density plots in Figl3](e) and (f) 
also support the proper description of the DFT+/7 for 
the structural transition of R-V02- As compared to the 
case of the DFT, the charge density in the DFT+/7 is 
seen to be more anisotropic, reflecting that more charges 
are accumulated in the bonding region along the c-axis. 
This phenomenon arises from the correlation-induced or- 
bital redistribution, whereby aig state is more occupied 
than other states. Accordingly, the system becomes 
more one dimensional-like along the c-axis, and so more 
susceptible to the Peierls transition. 

Now we have investigated the stability of M1-VO2 un- 
der pressure. In the DFT calculation for M1-VO2, we 
have obtained the insulating phase at the ambient pres- 
sure, even though the energy gap is negligibly small (~ 
0.03 eV). This result is different from existing DFT re- 
sults, in which metallic phases were obtained for Mi- 
V02.[l-[5| This difference is thought to come from the full 
relaxation of volume and ionic positions in our band cal- 
culations, which was not taken into account in the previ- 
ous calculations. Without the relaxation, we also get the 
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metallic phase for M1-VO2 in the DFT. In the DFT+[/, 
we have obtained the insulating phase with energy gap of 
~ 0.67 eV, which agrees well with the experimental gap 
of 0.6-0.7 eV. 0,113 

We have calculated the phonon dispersions of M1-VO2 
both in the DFT and the DFT+U. In both cases, there 
are no phonon softening instabilities, which indicate that 
M1-VO2 is a stable structure at the ambient pressure. 
Interesting feature is obtained in the phonon dispersion 
of M1-VO2 under pressure. Figure H] shows the phonon 
dispersion curve of M1-VO2 at 35.1 GPa in the DFT, 
which manifests the phonon softening instability at F. 
This softened mode implies that M1-VO2 would undergo 
a structural transition at some pressure P < 35.1 GPa. 
This result is indeed consistent with the experiments, [26- 
[isj in which the pressure-induced transition from insulat- 
ing M1-VO2 to a metallic phase of Mc-V02 was observed. 

The normal mode of the softened phonon at F is de- 
picted in the inset of Fig. HI which shows the rotation 
of oxygen (O) octahedron with slight displacements of V 
ions. The transformation from M1-VO2 to ME-VO2 is 
expected to arise from this rotation mode. In the experi- 
ment under pressure, the anisotropic compression was ob- 
served in Me-V02.[28] Lattice constants b and c become 
softened and hardened, respectively, while a decreases 
regularly. This anisotropy is reminiscent of the case of 
Cr02, which also has a pressure-induced structural tran- 
sition from rutile to CaCl2-type with the softened normal 
mode of the O octahedron rotation. Anisotropic com- 
pression is also present in Cr02 of CaCl2-type.[45|, |46[ 
This similarity suggests that ME-VO2 would have a mon- 
oclinic structure of CaCl2-type. 

In conclusion, based on the phonon dispersion stud- 
ies in the DFT and the DFT+/7, we have demonstrated 
that the driving mechanism of the MIT and the struc- 
tural transition in VO2 is the Mott-Peierls transition. 
The Coulomb correlation effect in R-YO2 plays an essen- 
tial role of assisting the Peierls transition. The softened 
phonon mode at q = R in the DFT+/7 describes prop- 
erly the structural transition from i?-V02 to M1-VO2. 
Further, we have found that M1-VO2 becomes unsta- 
ble under high pressure due to the phonon softening at 
F. This softened mode is expected to have the relation 
to the pressure-induced transition from M1-VO2 to M^- 
VO2, which is predicted to have a monoclinic structure 
of CaCl2-type. 
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